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Convolutional codes have been the centra/ part of most error-control systems in 
deep-space communication for many years. Almost all such applications, however, 
have used the restricted class of (n ,1), also known as {t rate 1/n, convolutional codes. 
The more general class of (n,k) convolutional codes contains many potentially useful 
codes , but their algebraic theory is difficult and has proved to be a stumbling block 
in the evolution of convolutional coding systems. In this article, the situation is 
improved by describing a set of practical algorithms for computing certain basic 
things about a convolutional code (among them the degree , the Forney indices, a 
minimal generator matrix, and a parity-check matrix), which are usually needed 
before a system using the code can be built. The approach is based on the classic 
Forney theory for convolutional codes, together with the extended Smith algorithm 
for polynomial matrices, which is introduced in this article. 


I. Introduction 

In his celebrated paper on the algebraic structure of 
convolutional codes, Forney [2] showed that by using the 
algebra of k x n polynomial matrices, in particular the 
invariant-factor theorem (aka the Smith Form), one can 
transform an arbitrary generator matrix for an ( n,k ) con- 
volutional code C into a noncat astrophic, basic, and ulti- 
mately minimal, generator matrix for C. He also showed 
how to find a polynomial inverse for a basic generator ma- 
trix for C, and a basic generator matrix for the dual code 
C 1 . In this article, efficient ways are discussed to do all 
these things. The main tool is an algorithm, called the 
extended Smith algorithm, which is used to find the in- 
variant factors of an arbitrary k x n matrix over an Eu- 


clidean domain, which bears the same relationship to the 
usual invariant factor algorithm as the extended Euclid’s 
algorithm bears to the usual Euclid’s algorithm. 

A brief review of Euclid’s algorithm (see e.g. [5, Section 
1.1] or [9, Chapter 2]) is helpful. The algorithm’s goal is to 
take a pair (a, b) of elements from a Euclidean domain R, 
and by repeated application of the division algorithm for R 
to find the greatest common divisor d of a and 6. The goal 
of the extended Euclidean algorithm (see e.g. [6, Section 
4.5.2] or [8, Section 8.4]) is to take the same pair, and not 
only find d, but also find elements s and t of R , such that 
sa+tb = d. The extended Euclidean algorithm has several 
important applications in coding theory. For example, it 
can be used to compute inverses in finite fields [9, Exam- 
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pie 4,2], decode BCH codes [8, Section 8.5], and to find 
finite impulse response (FIR) inverses for noncatastrophic 
(n, 1) convolutional generator matrices [1, Section 12,2], 
Similarly, the goal of the Smith algorithm (see e.g. [3, 
Section 6.2.4], [4, Section 6.3.3], or Smith’s nineteenth- 
century original article [10, Section 12.2]) is to take an 
arbitrary Jfexn matrix G (with k < n) over an Euclidean 
domain R , and by a sequence of elementary row and col- 
umn operations, to reduce G to a k x n diagonal matrix 
T = diag(7 i, . . . ,7*), whose diagonal entries are the in- 
variant factors of G, i.e., 7,* = A{/A t _i, where A i is the 
greatest common denominator of the i x i minors of G. 
(Here, A 0 = 1 is taken by convention.) Smith’s algorithm 
is reviewed in Section II. 

The goal of the extended Smith algorithm , which is in- 
troduced in this article, is to take the same input, and not 
only find T, but also to find a k xk unimodular matrix X, 
and annxn unimodular matrix Y , such that XGY = T. 
It is worthwhile to note that in the special case k = 1 
and n = 2, the (extended) Smith algorithm reduces to 
the (extended) Euclidean algorithm. The extended Smith 
algorithm is described in detail in Section III. 

Section IV describes how, given an arbitrary generator 
matrix G for an (n, J?) convolutional code, the extended 
Smith algorithm can be used to efficiently compute the 
things mentioned above (the degree, the Forney indices, a 
minimal generator matrix, a polynomial inverse, a parity- 
check matrix, etc., for the given code). Throughout this 
article, all results are illustrated with an example of a 2 x 4 
matrix of polynomials over the binary field GF(2), which 
is a generator matrix for a (4,2) convolutional code over 
GF(2). 


II. Smith's Algorithm 

In this section, a careful description is given of Smith’s 
algorithm (which is often called the invariant-factor algo- 
rithm), but a formal proof of its correctness is not given. 
For that, refer to [3, Section 6.2.4], [4, Section 6.3,3], or 
[11, Section 12.2]. 

Begin by recalling the definition of a general Euclidean 
domain [8, Chapter 2], It is an integral domain, i.e., a 
ring that satisfies the cancellation property, together with 
a “size” function |a| defined for every nonzero element a £ 
R. The size function must satisfy |a| < \ab\ if b ^ 0. 
Furthermore, if a and b are arbitrary, and 6 ^ 0, a can be 
“divided” by b in the sense that there exist elements q (the 
“quotient”) and r (the “remainder”) such that a — qb+r, 
where either r = 0 or else \r\ < |6|. For application to 
the study of convolutional codes, always take R to the 


ring of polynomials over a field F, where the “size” of a 
polynomial is its degree. However, there are many other 
Euclidean domains [8, Chapter 2], and Smith’s algorithm 
applies equally to all of them. 

The central part of Smith’s algorithm is the following 
subalgorithm E, which takes as input an arbitrary matrix 
A with entries from an Euclidean domain F, at least one 
of which is nonzero, and, via elementary row and column 
operations, transforms A into a matrix with the (1, 1) entry 
nonzero, every other entry in row 1 and column 1 zero, and 
every other entry in the matrix divisible by the (1,1) entry. 

El. Move the entry in A of least size to position (1,1). 

E2. If there is a nonzero entry in either row 1 or col- 
umn 1 that is not divisible by a n , use it to reduce 
the size of an, as follows: 

E2a. If aij, the entry in row 1 and column j, 
is not divisible by an, then by the di- 
vision algorithm there exist nonzero ele- 
ments q and r such that a^ = qa\i + r > 
with \r\ < |ai,i|. Thus, if q times col- 
umn 1 is subtracted from column j , the 
(l f j) entry will be changed to r, which 
has a smaller size than a \ r If this entry 
is moved to position (1, 1), the size of an 
will have been reduced. 

E2b. If an, the entry in row i and column 1, 
is not divisible by an, repeat the proce- 
dure outlined in step E2a, with rows and 
columns interchanged. 

E3. Reduce all the entries in row 1 and column 1 (ex- 
cept ai ( i itself) to zero, as follows: 

E3a. Since (by step E2) a xj for j > 2 is divisible 
by an, aij = qj^n- Thus, by subtract- 
ing qj times column 1 from column j , the 
(l,j) entry will be transformed to zero. 
Repeating this step for j = 2, . . . ,r, the 
first row will “zero out.” 

E3b. Repeat step E3a with rows and columns 
interchanged. 

E4. If there is a nonzero entry in A , say a^, which is 
not divisible by an, add column j to column 1. 
(This produces a nonzero entry in column 1 that 
is not divisible by a n .) Return to step E2. 

E5. Stop. The matrix A now has the desired property. 
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In the Smith algorithm itself, one applies subalgo- 
rithm E successively to the original matrix G, then to 
the matrix obtained from G by deleting row 1 and col- 
umn 1, then by deleting row 2 and column 2, etc., until 
either an all-zero matrix is encountered, or until all rows 
and columns have been processed. The result is a k x n 
diagonal matrix T = diag(7 x, . . . ,7*), whose diagonal en- 
tries are the invariant factors of G. The following example 
illustrates the Smith algorithm, when the underlying Eu- 
clidean domain is the ring of polynomials over the field 
GF(2), in which the “size” of an element is its degree. 


r -f 1 0 0 0 ^ 

5 "’VO 1 + D-fD 2 D + D 2 + D 3 1 + D 3 ) 

At this point, all of the remaining entries in row 2 are 
divisible by the (2, 2) entry: D+D 2 + D 3 = D (l+D+D 2 ) } 
and 1 -f D 3 = (1 + D) • (1 + D + D 2 ). Thus by adding D 
times column 2 to column 3, and then 1-f D times column 2 
to column 4, compute successively 

r _ (l 0 0 0 \ 

6 ~ \0 1-f D + D 2 0 1 + D 3 ) 


Example 1. Consider the following 2x4 matrix with 
entries that are polynomials over GF(2): 


f 1 0 0 (A 

^0 1 + D-fD 2 0 0 ) 


(\ 1 + D+D 2 1 + D 2 1 + D\ 
\D l + D+D 2 D 2 1 ) 


Since the lowest degree term already appears in the (1,1) 
position, skip step El. Since the (1,1) entry is 1, every 
nonzero entry in row 1 and column 1 is divisible by the 
(1,1) entry, so skip step E2 as well. Executing step E3: 
“zero out” the first row by adding 1 + D -f D 2 times col- 
umn 1 to column 2, 1 + D 2 times column 1 to column 3, 
and 1 + D times column 1 to column 4, thereby obtaining 
successively 


Gi 

G 2 

G 3 


a 

a 

a 


0 

1 + D 2 

1 + 0 

1 + D 3 

D 3 

1 

0 

0 


1 + D 3 

D + D 2 

+ D 3 

0 

0 


1 + D 3 

D + D 2 

+ D 3 


+ 

1 


D 


) 


0 

-fD + D 2 


To zero the entry in position (2, 1), add D times row 1 to 
row 2, thereby obtaining 


(10 0 0 

\0 1 + D 3 D+D 2 + D 3 1 + D + D 2 


This completes the operation of subalgorithm E on the 
original matrix G. Now apply subalgorithm E to the 1x3 
matrix obtained by deleting row 1 and column 1 from G 4. 
The lowest degree term in G4 appears in position (2, 4), so 
by interchanging columns 2 and 4, it is moved to position 
( 2 , 2 ): 


This completes the operation of the Smith algorithm on 
G. The matrix G 7 is the invariant-factor form for G; in 
particular, the invariant factors of G are 71 = 1 and 72 = 
1-f D + D 2 , ■ 

III. The Extended Smith Algorithm 

The previous section showed, in computational detail, 
how the Smith algorithm works. Beginning with the ma- 
trix Go = G, it produces a sequence of k x n matrices G,-> 
where G»+i is derived from G t by either an elementary 
row operation or an elementary column operation. This is 
represented algebraically as 

£*,-+1 = Ei+iGiFi+i ( 1 ) 

where and Fi+ \ are fcxfcandnxn elementary ma- 
trices, respectively. If G;+i is obtained from G, via a row 
operation, then \ is a nontrivial k x k elementary ma- 
trix, but F,*+i = I n . If Gi + i is obtained from G, via a 
column operation, then 1 is a nontrivial n x n elemen- 
tary matrix, but E i+ 1 = h. After a finite number N of 
steps, Gn = r. 

The extended Smith algorithm builds on the Smith al- 
gorithm. Whereas the Smith algorithm works only with 
the sequence Go, G\ , . . . , Gjv, the extended Smith algo- 
rithm also works with a sequence of unimodular k x Jc 
matrices Xo, . . . , A r J v , and a sequence of unimodular n x n 
matrices Yq, . . . , Y^. The sequences (X{) and (Y) are ini- 
tialized as Xo — Ik and Yo = / n , and updated via the rule 
[cf. Eq. (1)] 


X <+ 1 = Fj+iXj 


( 2 ) 
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( 3 ) 


If Eq. (4) is specialized with i = N, 


Yi+\ = YiFi + 1 

The following simple lemma is the key to the extended 
Smith algorithm. 

Lemma. 

XiGYi = Gi fori = 0,1,..., N (4) 


Proof. Since Xq = h and Yo = I n , Eq- (4) holds for 
i — 0. Assuming now that Eq. (4) holds for a given value 
of i, multiply both sides of Eq. (4) on the left by £’>+ 1 and 
on the right by I '\ + i 

Ei+yXiGYiFi+i = E i+1 GiF i+ i (5) 

By Eq. (1), the right side of Eq. (5) is equal to G,+i . Thus 

(E i+ iX,)G(Y;F i+1 ) = G i+ i (6) 

But E i+1 Xi = X i+1 by Eq. (2) and Y<F j+ 1 = Yi+i by 
Eq. (3), so Eq. (6) becomes simply Xi+\GYi+i = G{+ 1, 
which proves Eq. (4) by induction. I 


x n gy n = r (7) 

which is the desired extended Smith diagonalization of G. 
A convenient way to implement the extended Smith al- 
gorithm is to extend the original matrix G to dimensions 
(n + k) X (n + k) as follows: 



where in Eq. (8) 0 nx * is an n x k matrix of zeros. Then 
if the sequence of elementary row and column operations 
generated by the Smith algorithm applied to G is per- 
formed on the extended matrix G', after i iterations, the 
resulting matrix GJ has the form 



Thus, after N steps, the matrices A' ; y and V'y in Eq. (7) 
appear as the upper-right and lower-left blocks of G' N , re- 
spectively. All this is illustrated by extending the example 
from Section II. 


Example 2. Begin as in Example 1, with the following 2x4 matrix over GF(2)[D]: 


^ ~ _ ( 1 1 + D + D 2 1 + Z> 2 1 + D\ 
G = Go ~\D l + D + D 2 D 2 1 ) 


Then the corresponding matrix G' [cf. Eq. (8)] is 


°\ 

1 

0 

0 

0 

0 / 


G' = G' = 


1 

0 

0 

Vo 


1 + D + D 2 
1 + D+ D 2 
0 
1 
0 
0 


1 + D 2 

D 2 

0 

0 

1 

0 


1 + D 
1 
0 
0 
0 
1 


1 

0 

0 

0 

0 

0 


In Example 1, Gi, G2, and G3 are obtained from Go by successively adding 1 + D + D 2 times column 1 to column 2, 
1 + D 2 times column 1 to column 3, and 1 + £> times column 1 to column 4. If the same sequence of operations is 
performed to Gq, one obtains successively 
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/I 

0 

1 + D 2 1 + D 

1 





D 

1 + D 3 

D 2 

1 

0 1 





1 

1 + D + D 2 

0 

0 

0 0 





0 

1 

0 

0 

0 0 





0 

0 

1 

0 

0 0 





V 0 

0 

0 

1 

0 0/ 





1 

0 

0 


1 + D 

1 




D 

1 + D 3 

D + D 2 + D 3 

1 

0 

1 



1 

1 + D + D 2 

1 + D 2 


0 

0 

0 



0 

1 

0 


0 

0 

0 



0 

0 

1 


0 

0 

0 



V 0 

0 

0 


1 

0 

0/ 



1 

0 

0 


0 



1 

°\ 

D 

1 + D 3 

D + D 2 + D 3 

1 + D + D 2 0 

1 

1 

1 + D + D 2 

1 + D 2 


1 + D 

0 

0 

0 

1 

0 


0 


0 

0 

0 

0 

1 


0 


0 

0 

0 

0 

0 


1 


0 

0/ 


Next, adding D times row 1 to row 2, as was done to obtain from G3, obtain 



1 

0 

0 

0 

1 

°\ 


0 

1 + D 3 

D + D 2 + D 3 

1 + D + D 2 

D 

1 


1 

0 

1 + D + D 2 
1 

1 + D 2 
0 

1 + D 
0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 


\0 

0 

0 

1 

0 

0 / 

Interchanging columns 2 and 4, obtain 







/I 

0 

0 

0 

1 

°\ 


0 

1 + D + D 2 

D + D 2 + D 3 

1 + D 3 

D 

1 

G' 5 = 

1 

0 

1 + D 
0 

1 + D 2 
0 

1 + D + D 2 
1 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 


\0 

1 

0 

0 

0 

0 / 


Finally, adding D times column 2 to column 3, and 1 + D times column 2 to column 4, one computes successively 


/I 

0 

0 

0 

1 

°\ 

0 

1 + D + D 2 

0 

1 + D 3 

D 

1 

1 

1 + D 

1 + D 

1 + D + D 2 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

\0 

1 

D 

0 

0 

0/ 


and 
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/I 

0 

0 

0 

1 

°\ 


0 

1 + D + D 2 

0 

0 

D 

1 

G' 7 = 

1 

0 

l + D 

0 

l + D 
0 

D 

1 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 



1 

D 

l + D 

0 

0/ 


Thus the extended Smith decomposition of the original matrix G is given by [cf. Eq. (7)] 

1 l + D l + D D \ 

0 0 0 1 | _ f l 0 0 0 \ 

0 0 1 0 ” VO 1 + D + D 2 0 0 ) 

0 1 D 1 + Dl 

In Section V, the decomposition of Example 2 will be used in the analysis of the (4,2) convolutional code generated by 
G, as an illustration of the general results to be expounded upon in Section IV. 



IV. Application to the Analysis of 
Convolutional Codes 

In this section, given an arbitrary generator matrix G 
for an (n,fc) convolutional code G, one can efficiently find, 
among other things; a basic generator matrix, say Gbasic 
for G; a polynomial inverse for Gbasic; a minimal generator 
matrix for G; and a basic generator matrix for the dual 
code C -1 ". The central tool needed to do all this is the 
extended Smith algorithm introduced in Section III. 

Assume then that G is a k x n generator matrix for 
an (n,fc) convolutional code over the field F, i.e., a k x n 
matrix of rank k over the field F(D) of rational functions 
over F . By multiplying the ith row of G by the least com- 
mon multiple of the denominators in that row, one easily 
obtains an equivalent generator matrix, all of whose en- 
tries are polynomials over F, i.e., elements of F[D], Since 
F[D] is an Euclidean domain, one may apply the extended 
Smith algorithm described in Section III, thereby obtain- 
ing a decomposition of Eq. (7). In what follows, the matri- 
ces Xn and Y Xr produced by the extended Smith algorithm 
will be denoted simply by X and Y. 

The matrices X, Y, and V contain much valuable in- 
formation about the code C and the generator matrix G. 
To extract this information, however, first define several 
useful pieces of these matrices, called r*, T*, K ) and H 
(in what follows, r = n — k) 

Th = leftmost k columns of T = diag(7 i, . . . ,jk) 

(dimensions: k x k) (10) 


ft = 7 * • r* 1 = diag(7t/7i,. ..,7t/7i) 

(dimensions: k x k) (11) 

K = leftmost k columns of Y 

(dimensions: n x k) (12) 

H = rightmost r columns of Y 

(dimensions: n x r) (13) 

The following theorem describes the useful outputs of 
the extended Smith algorithm when applied to G. (For 
terminology not fully explained here, refer to Forney [2].) 

Theorem 1. With the matrices T*, f*, 7i, and H 
defined as in Eqs. (10)— (13), one has the following: 

(a) A basic generator matrix for G : Gbasic = T^ 1 A r G. 
(That is, Gbasic is obtained by dividing the fth row of 
XG by the invariant factor 7 », for i — 1 , . . . , k.) 

(b) A polynomial inverse for Gbasic' 

(c) A polynomial pseudo-inverse for G, with factor 7*: 
ATfcX. (In particular, if G is already basic, i.e., if 
T* = /*, then ATX is a polynomial inverse for G.) 

(d) A basic generator matrix for C L , i.e. } a parity-check 
matrix for C: H T , 


27 



Proof. From Eq. (7) with X = X N and Y -Y N) one 
has 

XG — TB (14) 

where B = Y -1 . Now B is an nxn unimodular polynomial 
matrix. Let V be the k x n matrix consisting of the first 
k rows of B , and L be the r x n matrix consisting of the 
last r rows of B. Then, since T = ( IT 0 * X r ) where Ofc Xr 
is a A: x r all-zeros matrix, it follows that 

TB = (T k <W)(5Q=I»17 (15) 

so that, combining Eqs. (14) and (15) 

(16) 

From Eq. (16), U is row equivalent to G. Furthermore, 
since U = {h Ofcxr ) B, and B is unimodular, it follows 
that the invariant factors of U are all equal to 1 , and so U is 
a basic generator matrix for C. This proves Theorem 1(a). 

To prove Theorem 1(b), use the fact that Y and B are 
inverse matrices, so that 

-= 00 <* *>-(« ») 

='•=(»!:, °r) < i7 > 

It follows from Eq. (17) that UK = , so that K is a 

polynomial inverse for the basic generator matrix U> which 
proves Theorem 1(b). 

To prove Theorem 1(c), suppose that J is an n x k 
polynomial pseudo- inverse for G with (polynomial) factor 
^(D), i.e., that GJ = tpl k . Then 

XGYY~ 1 JX~ l = XGJX- 1 = Xtyl^X- 1 = iph (18) 

But from Eq. (7), XGY = T = (I* 0 txr ). Thus, if 
the matrix Y~ l JX~ 1 is denoted by J' , J' is a polyno- 
mial matrix (since X and Y are unimodular) and Eq. (18) 
becomes 

(r* o )j' = i>i k (19) 


If now J' = ^ where Jq is k x k and J[ is r x Ar, 

it follows from Eq. (19) that Jq = ^r^ 1 . But since J ! 
is a polynomial matrix, and T * = diag( 7 f , . . . ), it 

follows that ^ must be a multiple of 7 *. Conversely, if 
tp = 7 *, so that /q = TibT ^ 1 = T k and J[ — O rx Jfc> lh en 

J 1 = f J k \ Thus, J = YJ f X = ATjtX is a polynomial 

VW/ 

pseudo-inverse with factor 7 ^, which proves Theorem 1(c). 

To prove Theorem 1(d), one has to show that the set 
of codewords in the code generated by G, i.e., the set of 
n-dimensional E(D)-vectors y of the form y — xG , is iden- 
tical to the set of vectors y such that yH — 0. Thus, let 
Yi = {y : y - xG } and Y 2 = {y : yH ~ 0). It will be 
shown that Yi = Y 2 . 

Since, as in Theorem 1(a), the matrix {7, defined to be 
the first k rows of B = Y" 1 , is (a basic generator matrix) 
equivalent to G, it follows that Y\ = {y : y — xU) . Thus, 
if y G Yi, y = xU for some x and so yH — ( xU)H = 
x{U H) = 0, since by Eq. (17), UH = O^xr • Thus, Y\ C Y 2 . 
On the other hand, if y is arbitrary, then since B~ l = Y, 

y = yI n = y(YB) = y((K 

= y(I<U + HL) = (yK)U + ( yH)L (20) 

Now suppose y E Y 2 , i.e., yi/ = 0. Then from Eq. (20) 
y = ( yK)U — xt/, where x = yA", and so y G Y x . Thus, 
Yi C Y 2) which completes the proof of Theorem 1(d). I 

In the next section these results will be briefly illus- 
trated with a simple example. 


V. An Example Convolutional Code 

To illustrate the results of this article, consider the fol- 
lowing generator matrix G for a binary (4, 2) convolutional 
code: 

n _( \ 1 + D + D 2 1 + D 2 1 + D) 

° “ \D 14 - D + D 2 D 2 1 J 

In Section III the extended invariant-factor decomposition 
of G was found to be XGY = T, where 
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1 0 0 0 

0 1 + D + D 2 0 0 


1 + D 

1 + D 

D 

0 

0 

1 

0 

1 

0 

1 

D 

1 + D 


Thus, the matrices defined in Eqs. (10)— (13) are as follows: 


r *~(o 1 + D + D 2 ) 


1 + D + D 2 O' 
0 1 


I<T k X = 


1 1 + D 

0 0 

0 0 

D 1 


(d) A (basic) generator matrix for C A 


T _ (1 + D 0 1 D 


H 1 = 


D 1 0 1 + D 


Finally, note that in this case neither of the basic gener- 
ator matrices Gbasic or ls To minimize them, 

follow the simple algorithm originally described in [2], or 
perhaps more lucidly in Kailath [4, Section 6.3.2] (where 
minimal matrices are, however, called row-reduced). The 
idea is to use elementary row operations on G to reduce 
the degree of the highest degree terms in some row, as 
long as this is possible. For example, to minimize Gbasic? 
multiply the second row by D and add it to the first row, 
obtaining 


1 0 0 0 \ 

1 + D 0 0 \) 

1 + D 0 1 D 

D 1 0 1 + D 


Using the prescriptions in Theorem 1, one quickly ob- 
tains the following: 

(a) A basic generator matrix for C : 

Gbasic = rj'XG 


1 l + D + D 2 1 + D 2 1 + D 
0 1 + D D 1 


(b) A polynomial inverse for Gbasic* 


1 1 + D 

0 0 

0 0 

0 1 


11 11 

0 1 + D D 1 


Since elementary row operations cannot reduce the row de- 
grees further, the matrix G ! is a minimal generator matrix 
for the code G. (Indeed, at this point, one can recognize C 
as the df ree = 4 partial-unit-memory (4,2) code first found 
by Lauer [7, Table 1].) To find a polynomial inverse for 
G r , note that G' = TG, where T is a unimodular k x k 
polynomial matrix, and so K f = KT“ l is a polynomial 
inverse for G'. In this example, simply multiply the first 
column of K by D and subtract it from the second column , 
thereby obtaining the following polynomial inverse for the 
minimal generator matrix G': 


To minimize H T , simply add row 2 to row 1, thereby 
obtaining the following minimal generator matrix for the 
dual code: 


111 1 

D 1 0 1+D 


(c) A polynomial pseudo-inverse for G, with factor 72 
l + D + D 2 : 


In this case, the dual code is isomorphic to the original 
code. 
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The Forney indices of an (n,fc) convolutional code are 
the degrees of the rows of a minimal generator matrix, and 
the degree of the code is the sum of the Forney indices. 
In this example, the degree of both C and C **■ is 1, and 


the Forney index of both codes is (0,1). It is a general 
theorem that deg C = deg C L (see [2, Theorem 7]), but the 
equality of the Forney indices in this case is more or less 
accidental. • 
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